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ABSTRACT 

In this work we extend and improve the results done in a 
previous work on simulating Spiking Neural P systems (SNP 
systems in short) with delays using SNP systems without 
delays. We simulate the former with the latter over sequen- 
tial, iteration, join, and split routing. Our results provide 
constructions so that both systems halt at exactly the same 
time, start with only one spike, and produce the same num- 
ber of spikes to the environment after halting. 



Keywords: Membrane Computing, Spiking Neural P sys- 
tems, delays, simulation, routing 

1. INTRODUCTION 

SNP systems, first presented in [7] and with some recent 
results in [12], [13] and [10] (among others), are computing 
devices inspired by how biological neurons represent infor- 
mation: using electro-chemical signals called spikes. Since 
spikes are indistinct, information is taken not from the spikes 
themselves, but from their multiplicity or time of arrival. 
One motivation for SNP systems (as is the case in the area 
of Membrane Computing 11 in general) is to abstract ideas 
from biology for computational use. For SNP systems in 
particular, the neuron from our brains is the motivation. 
It can be argued that the human brain is one (if not cur- 
rently) the most complicated and powerful "supercomputer" 
known to us at the moment. The brain performs complex 
computations from billions of interconnected neurons while 
consuming only around 10 to 20 Watts of energy [8], and 
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it is small enogh to fit in our skulls. It is therefore desir- 
able to work with as little quantity of "energy" as possible, 
and we can think of the spike in SNP systems as being such 
quantity. 

SNP systems are Turing complete devices [5] and have 
been used as (among others) transducers [5], generating vec- 
tors of numbers Q~], as well solving hard problems 9 . Spik- 
ing rules (rules that produce spikes) are usually of two types: 
with delays and without delays. If an SNP system has at 
least one rule with a delay, we refer to it as an SNP system 
with delay labeled as II, otherwise it is known as an SNP 
system without delay labeled as II. 

In this work we extend the work presented in [3], with the 
goal of simulating a II that peforms sequential, iteration, 
join, and split routing with a II that performs the same 
routings. By routing we mean the transfer or movement of 
spikes from one neuron to another. By simulation in [3] it 
is meant that the following two requirements are satisfied: 

Ri : Halting time of n coincides with the halting time of n, 
or is offset either by a fixed timestep or by a function 
of the delays in n, 

R2 : number of spikes in the final configuration of n is the 
same number in n, or is offset by a function of the 
delays in n. 

In [3], the construction of n from n is such that the initial 
spikes of n is a function of the delay (or delays) in n. In 
particular, the initial spikes of n are multiples of the delays 
in n. Aside from the increased initial spike number in n, 
the exponents in the regular expressions and the number 
of consumed spikes of certain spiking rules in n are also 
multiples of all the delays in a given routing of n. 

We improve the work done in [3] by providing alternative 
constructions in this work. Our specific contributions are as 
follows: 



• we construct a n that simulates a n that performs 
sequential, iteration, join, and split routing, 

• both n and n start with only one spike each in the 
initial neuron, 

• halting time of n and n coincide i.e. there are no 
offsets, 



number of spikes sent to the environment after halting 
are equal for II and II. 

our construction allows split routing even if the delays 
of the output neurons are not equal. 



The trade-off is that for every delay d in II, we add d neurons 
in IT. If the initial neuron of II has a delay, following our con- 
struction means we simply modify II and II such that their 
new halting time involves one additional time step. The 
succeeding sections are as follows: Section [2] provides pre- 
liminaries and assumptions for our work. Section [3] presents 
our main results. We end with our final remarks and direc- 
tions for future work in Section [4] 

2. PRELIMINARIES 

It is assumed that the readers are familiar with the basics 
of Membrane Computing (a good introduction is [1 1 J with 
recent results and information in the P systems webpage 
at http://ppage.psystems.eu/ and a recent handbook in 
113] ) and formal language theory. We only briefly mention 
notions and notations which will be useful throughout the 
paper. 

Let V be an alphabet, V* is the free monoid over V with 
respect to concatenation and the identity element A (the 
empty string). The set of all non-empty strings over V is 
denoted as V + so V + = V* — {A}. We call V a single- 
ton if V = {a} and simply write a* and a + instead of {a*} 
and {a + }. The length of a string w 6 V is denoted by 
\w\. If a is a symbol in V, a — A. A language L C V* is 
regular if there is a regular expression E over V such that 
L(E) = L. A regular expression over an alphabet V is con- 
structed starting from A and the symbols of V using the 
operations union, concatenation, and +, using parentheses 
when necessary to specify the order of operations. Specifi- 
cally, (i) A and each a € V are regular expressions, (ii) if Ei 
and E2 are regular expressions over V then (E1UE2), E\E 2 , 
and E+ are regular expressions over V, and (Hi) nothing 
else is a regular expression over V. With each expression 
E we associate a language L(E) defined in the following 
way: (i) L(\) = {A} and L(a) = {a} for all a £ V, (ii) 
L(Ei U E 2 ) = L(Ei) U L(E 2 ), L(EiE%) = L(E t )L(E 2 ), and 
L(E^) = L(Ei) + , for all regular expressions E%, E2 over V . 
Unnecessary parentheses are omitted when writing regular 
expressions, and E + U {A} is written as E* . Next we have 
the definition for an SNP system. 



Definition 1 (SNP system). An SNP system, of a fi- 
nite degree m > 1 is a construct of the form 



II = (0,ai,..., cr m , syn,out), 



where: 



(b) Ri is a finite set of rules of the general form 



E/a ->• a b ;d 

where E is a regular expression over O, c > 1, if 
b > then d > and c > b, else if b = then 
d = 0. 

3. syn C {1, 2, . . . , m} x {1, 2, . . . , m}, (i, i) ^ syn for 
1 < i < m, are synapses between neurons. 

4- out € {1, 2, . . . , m} is the index of the output neuron. 



A spiking rule is where b > 1. A forgetting rule is a rule 
where b = is written as E/a° — > A. If L(E) = {a c } 
then spiking and forgetting rules are simply written as a c — > 
a b and a c — > A, respectively. Applications of rules are as 
follows: if neuron ai contains k spikes, a k £ L(E) and k > c, 
then the rule E/a c — > a b G Ri is enabled and the rule can 
be fired or applied. If b > 1, the application of this rule 
removes c spikes from ai , so that only k — c spikes remain in 
Oi. The neuron sends b number of spikes to every en, such 
that (i,j) € syn. The output neuron has a synapse not 
directed to any other neuron, only to the environment. The 
neuron o\ is referred to as the initial neuron. 

If a spiking rule (forgetting rules cannot have delays) has 
d — 0, the b number of spikes are sent immediately i.e. in 
the same time step as the application of the rule. If d > 1 
and the spiking rule was applied at time t, then the spikes 
are sent at time t + d. From time t to t + d — 1 the neuron 
is said to be closed (inspired by the refractory period of the 
neuron in biology) and cannot receive spikes. Any spikes 
sent to the neuron when the neuron is closed are lost or 
removed from the system. At time t + d the neuron becomes 
open and can then receive spikes again. The neuron can 
then apply another rule at time t + d + 1. If 6 = then no 
spikes are produced. SNP systems assume a global clock, so 
the application of rules and the sending of spikes by neurons 
are all synchronized. 

A configuration of the system at time k is denoted as Ck = 
(ni/ti, . . ., n m /t m ,n e ), where each element of the vector 
(except for n e , denoting the spikes in the environment) is 
the configuration of a neuron Gi, with nt spikes and is open 
after ti > steps. An initial configuration Co is therefore 
(ni/0, . . . , n m /0, 0) since no rules whether with or without 
delay, have yet been applied and the environment is ini- 
tially empty. A computation is a sequence of transitions 
from an initial configuration. A computation may halt (no 
more rules can be applied for a given configuration) or not. 
If an SNP system does halt, all neurons should be open. 
Computation result in this work is obtained by checking the 
number of spikes in the environment once the system halts. 



1.0 = {a} is the singleton alphabet (a is called spike). 

2. (Ti, ... , Om are neurons of the form Oi = (n;, Ri), 1 < 
i < ra, where: 

(a) Hi > is an integer representing the number of 
spikes in o~i 



—a x s 

/a — > a ) Ja + /a 



O; 2 j J o / a 




Figure 1: SNP system with delay IIo. 



As an example, let us have an SNP system shown in Figure 
[T] formally defined as follows: Ilo = (O, o~\, 02, 03, syn, out) 
where o\ — (l,a + /a — > a), o 2 = (0,a + /a — > a; 2), 03 = 
(0, a + /a — > a), syn — {(1,2), (2,3)}, initial neuron is <7i, 
and out — 3. Only neuron 01 has one spike initially and 
only 02 has a rule with a delay d = 2. We have Co = 
(1/0,0/0,0/0,0). At the next step, 01 can use its rule (it 
has at least one spike) and consumes one spike and sends one 
spike immediately to 02 so we have Ci = (0/0, 1/0,0/0,0). 
At step 2, (72 consumes its spike and closes for 2 time steps, 
so C* 2 = (0/0, 0/2, 0/0, 0). At step 3 we have C 3 = ( 0/0, 0/1, 
0/0, 0). At time step 4, 02 opens and sends one spike to 03, 
so Cz — (0/0, 0/0, 1/0, 0). Finally, at time step 5 the output 
neuron sends one spike to the environment, llo halts and we 
have C 5 = (0/0,0/0,0/0,1). 




Figure 2: Routing constructs (from left to right): 
sequential, split, and join. 

SNP systems where each neuron has exactly one rule are 
called simple, while the systems that have the same set of 
rules are called homogeneous [14] . In this work, if SNP sys- 
tems only have rules of the restricted form (a k ) + /a k —¥ a 
where k is a non-negative integer, we refer to them as semi- 
homogeneous. We only consider SNP systems II and II that 
are simple and semi-homogeneous, where their initial con- 
figurations have one spike in the initial neuron only, and no 
spike in every other neuron (as in Figure [TJ in this work. We 
make no restrictions on the values of the delays in a (rule 
of a) neuron. The objective is to route or move the single 
spike in the initial neuron through the system, towards the 
output neuron, and eventually to the environment. Spikes 
are routed via paths, where a path consists of at least two 
neurons &i,Oj such that (i,j) £ syn. Using paths, we can 
have four basic routing constructs (referring to Figure [2]): 

1. sequential where, given at least two neurons o\, 02 
such that o 2 spikes only after o\ spikes and there is a 
path from o\ to o 2 , 

2. iteration, where at least two neurons spike multiple 
(possibly an infinite) number of times and a loop is 
formed e.g. adding a synapse (2, 1) which creates a 
loop between o\ and o 2 , 

3. spltt, where a spike from 0-3 is sent to at least two 
output neurons 04 and (75 and (3, 4), (3, 5) £ syn, 

4. join, where spikes from at least two input neurons 
06,07 are sent to a neuron erg, where (6, 8), (7, 8) G 
syn, so that as produces a spike only after accumulat- 
ing spikes from or and o$. 



Notice that iteration routing can be formed by combining 
the three other constructs. Also notice that if there exists 
a sequential path from Oi (with delay di) to Oj (with delay 
(fe) so that di < d2 and the number of spikes of the initial 
neuron o\ in Co is n% > 1, it is possible for some spikes to 
be lost. The reason is that it is possible for Oj to still be 
closed when spikes from Oi arrive. We avoid lost spikes by 
considering SNP systems where the initial neuron has only 
one spike. We say in this work that a II simulates a II if two 
requirements are satisfied: 

R'i : halting time of II spikes is the same halting time of II, 

R2 : number of spikes in the environment of II when II halts 
is equal to the number of spikes in the environment of 
II when II halts. 

3. MAIN RESULTS 

We begin presenting our results with a fundamental idea on 
sequential routing. 

Lemma 1 (Sequential routing). Given an SNP sys- 
tem with delay II performing sequential routing, there exists 
an SNP system without delay II peforming sequential routing 
that simulates II. 

Proof. We refer to Figure[3]for illustrations. Let 

II = (0,cm,cri2, {(11, 12)}, 12) with (7i = {l,a+/a a) 
and o"2 = (0, a + /a — > a; d) 

we then let 

II = (0, 021, 022-i, (723, syn, 23) where 1 < i < d, syn — 
{(21,22-1), ...,(21,22-(d- 1)), (22-1,22-d), . . . , (22-(d - 
1), 22-d), (22-d, 23)}, a 2 i = (1, a+/a -)■ a), 022-i = (0, a+/a -> 
a), (723 = (0,(a d - 1 ) + /a d - 1 ->a). 

The additional d neurons, immediately after initial neuron 
(72i in LI, are used to multiply the single spike from 021. 
The additional neurons then send one spike each to 022-d ■ 
Neuron 022-d accumulates d — 1 spikes, and consumes these, 
one spike at a time and sending one spike every time to 
(723. This consumption of one spike every time step creates 
a delay of d — 1 time steps. Due to the regular expression 
of the rule in 023 , the neuron will have to accumulate d — 1 
spikes before the rule is used. Once 023 accumulates d — 1 
spikes, it immediately sends one spike to the environment. 
This spiking and halting occurs at time t + d + 1 for both II 
and n (satisfying Ri and R2) if we let t be the time when 
<7n and (721 spike. 

We can repeatedly apply the previous construction if there 
exist more than one neuron with a (rule having a) delay 
in a sequential path as seen in Figure [3] It can be easily 
shown that if there exists Oi without delay in a sequential 
path between o\\ and (712, the time to halt for both IIi and 
IIi still coincide. In particular, every additional Oi having a 
rule without a delay adds one time step to the halting time 
of both IT and II. Both R^ and R 2 are still satisfied. □ 




Figure 3: Sequential routing: IT (top) with delay d, 
and IIi (bottom) simulating IIi. 

A sample computation of IIi and IIi is shown in Table [1] 
For sample computations of II2 and IT 2 we refer to Tableful 
From Lemma [T] we have the following observation. 

Observation 1. If IT has more than one neuron with a 

?n 

delay in a rule, the total additional neurons in II is di 

i=i 

where di is the delay of the rule in a . 



Steps 


Hi 


III 


to 


(1/0,0/0,0) 


(1,0,0,0,0,0) 


h 


(0/0, 1/0, 0) 


(0,1,1,0,0,0) 


ta 


(0/0,0/3,0) 


(0,0,0,2,0,0) 


ts 


0,0 2,0 


(0,0,0,0,1,0) 


U 


(0/0,0/1,0) 


(0,0,0,0,2,0) 


ts 


(0/0,0/0,1) 


(0,0,0,0,0,1) 



Table 1: Sample computations of IIi and IIi, d — 3. 



Steps 


n 2 


n 2 


to 


(1/0, 0/0, 0/0, 0) 


(1, 0, 0, 0, 0, 0, 0, 0, 0) 


ti 


(0/0, 1/0, 0/0, 0) 


(0, 1, 0, 0, 0, 0, 0, 0, 0) 


ts 


(0/0, 0/2, 0/0, 0) 


(0, 0, 1, 0, 0, 0, 0, 0, 0) 


ts 


(0/0, 0/1, 0/0, 0) 


(0, 0, 0, 1, 0, 0, 0, 0, 0) 


U 


(0/0, 0/0, 1/0, 0) 


(0, 0, 0, 0, 1, 1, 0, 0, 0) 


t 5 


(0/0, 0/0, 0/3, 0) 


(0, 0, 0, 0, 0, 0, 2, 0, 0) 


te 


(0/0, 0/0, 0/2, 0) 


(0, 0, 0, 0, 0, 0, 1, 1, 0) 


h 


(0/0, 0/0, 0/1, 0) 


(0, 0, 0, 0, 0, 0, 0, 2, 0) 


t 8 


(0/0, 0/0, 0/0, 1) 


(0, 0, 0, 0, 0, 0, 0, 0, 1) 



Table 2: Sample computations of II2 and IT2, di = 2, 
d 2 = 3. 



Lemma 2 (Iteration routing). Given an SNP sys- 
tem with delay II performing iteration routing, there exists 
an SNP system without delay IT peforming iteration routing 
that simulates IT. 

Proof. We refer to Figure [S] for illustrations. Let 

n = (0,0-n, o-ia, {(11, 12), (12,11) },12) where an = (1, 
a + /a —} a; d, 02 = (0, a + /a — > a) 

we then let 




25 \ 24-d 2 



Figure 4: Sequential routing with multiple delays: 
IT2 (top) with delays di and d2, and IT2 (bottom) 
simulating IT2. 

IT = (0,a 2 i,a 2 2-i, 023, syn, 23) where 021 = (1, a + /a — > a), 
0-22-i = (0, a+/a -+ a) for 1 < i < d, a 23 = (0, {a d - 1 ) + /a d - 1 - 
a) and syn = {(21, 22-1), . . . , (21, 22- (d - 1)), (22-1, 22-d), 
. . . , (22- (d - 1), 22-d), (22-d, 23), (23, 21)}, out = 23 

The construction of LI uses the construction idea in Lemma 
[T]i.e. the neuron with a delay d in IT is replaced with d addi- 
tional neurons in IT. In Figure [S] an infinite loop is created: 
a spike starts at an and it uses its rule at time t so that 
the spike is sent to ai2 at time t + d, then ai2 immediately 
sends a spike back to an (and the environment) at time 
t+d+l, and so on and so forth. Similarly, a2i sends a spike 
to neurons C22-1 to a 2 2-(d-i) at time t. At time t+1, 022-d 
accumulates d — 1 spikes from the d — 1 neurons from the 
previous time step. The spikes in a22-d are consumed and 
then sent one at a time to a23 . At time t + d, 023 accumu- 
lates d — 1 spikes so that it sends one spike back to 021 and 
at the environment at time t + d+l, coinciding with the 
time of spiking of ai2. Thus, and R 2 are satisfied. □ 

Lemma[2]for iteration routing makes use of the construction 
used in Lemma [T] for sequential routing. This construction 
will again be used for the join and split routings as follows. 
From Lemma [2] we have the following observation. 

Observation 2. If the initial neuron of IT has a delay 
and its halting time is t + d, we add a new initial neuron 
ay in IT with (1', 1) £ syn so that IT halts at time t + d+l. 
We then add a new initial neuron similarly to IT and modify 
its syn (following Lemma\T\ construction) so that IT halts at 
t + d+l instead, simulating IT. 

Although the premise of Observation [2] is different from our 
assumption in Section [2] that the initial neuron has no de- 
lay, the observation provides a solution on how to approach 
such a premise. For example, if an has a delay instead 
of a23 in IT3, we add a new initial neuron to a iv and a 



new synapse (11', 11) € syn in II. For II3, we modify it as 
follows: add a new initial neuron o 2 i' an d (T21 is replaced 
with d neurons (instead of 0-22). The synapse set of II is 
changed to {(21', 21-1), . . . , (21', 21-(d-l)), (21-1, 21-d), . . ., 
(21-(d- 1), 21-d), (21-d, 22), (22, 21-1), . . . , (22, 21-(d- 1))}, 
we remove 023 and have (J22 as the output neuron instead. 
Both II3 and 1T3 halt at the same time at t + d + 2. 




Figure 5: Iteration routing: II3 (top) has a delay d 
simulated by II3 (bottom). 

Lemma 3 (Join routing). Given an SNP system with 
delay II performing join routing, there exists an SNP system 
without delay II peforming join routing that simulates II. 

Proof. We refer to Figure [6] for illustrations. Let 

II = (0,0-11,0-12,0-13, {(11, 13), (12, 13)}, 13) where an = 
CT12 = (1, a + /a — ► a), 1T13 = (0, (a 2 ) + /a 2 —¥ a; d) 

we then let 

II = (O, 021, 022, 023-1, C24, syn, 24) where 1 < i < d, an = 
022 = (l,a + /a ->■ a), 023-1 = ... = 2 3-(d-i) = (0, (a 2 ) + /a 2 -> 
a 2 ), 023-d = (0,(a 2 ) + /a 2 -> a), syn = {(21,23-1), 
(21,23-(d-l)), (22,23-1), (22, 23-(d- 1)), (23-1, 23-d), 
. . . , (23-(d - 1), 23-d), (23-d, 24)}. 

For II and II we have as initial neurons an, 012 and 0-21 , 0-22 
respectively. Using the construction in Lemma [T] II has d 
additional neurons corresponding to 013 in II. Let time t be 
the time when the initial neurons spike. At time t, neurons 
CT23-1 to 023-(d_i) have two spikes each, so that in total, II 
at this time has 2(d — 1) spikes. At the next time step t + 1, 
the d — 1 neurons send two spikes each to 0-23-d so that 023-d 
accumulates 2(d — 1) spikes. Since 023-d consumes two spikes 
every time and it has 2(d — 1) spikes, 0-23-d will take d — 1 
time steps to consume all of its 2(d — 1) spikes. Every time 
023- d spikes, it sends only one spike to 024. At time t + d, 024 
has accumulated d — 1 spikes from 023-d so that 024 sends 
one spike to the environment and halts at t + d + 1. This 
time step coincides with the halting time of 013 in II, sending 
one spike to the environment. We therefore satisfy R^ and 
R 2 . □ 

A sample computation of II4 and II4 is shown in Table [3] 




Figure 6: Join routing : II4 (top) has delay d simu- 
lated by II4 (bottom). 



Steps 


n 4 


n 4 


to 


(1/0,1/0,0/0,0) 


(1,1,0,0,0,0,0) 


h 


(0/0, 0/0, 2/0,0) 


(0,0,2,2,0,0,0) 


ti 


(0/0, 0/0,0/3,0) 


(0,0,0,0,4,0,0) 


t 3 


(0/0, 0/0,0/2,0) 


(0,0,0,0,2,1,0) 


u 


(0/0, 0/0,0/1,0) 


(0,0,0,0,0,2,0) 




(0/0, 0/0,0/0,1) 


(0,0,0,0,0,0,1) 



Table 3: Sample computations of II4 and II4, d = 3. 



Lemma 4 (Split routing). Given an SNP system with 
delay II performing split routing, there exists an SNP system 
without delay II peforming split routing that simulates IT. 

Proof. We refer back to the split routing in Figure [5] 
Notice that a split routing can be thought of as two separate 
paths, either from 0-3 to 04 or 03 to 0-5 . We let t be the time 
that 03 spikes and modify the split routing in Figure [2] as 
follows and let it be 

II = (0, 0-3, 04, 0-5, a , syn, o) where 03 = (l,a + /a — > a), 
05 = (0, a + /a — > a), 04 = (0,a + /a — > a;d), syn = {(3,4), 
(3,5), (4,o), (5,o)}. 

We arbitrarily chose 04 to have a delay instead of 05 in this 
case. Next we let 

II = (O, 03/ , 04/.j, o 5 / , a a , syn,o) where 1 < i < d, 03/ = 

(l,a + /a ->■ a), a A '-\ = ... = oV-fd-i) = °V = (0,a+/a 

a), <T4>-d = (0,(a d - 1 ) + /a d - 1 -». a), syn = {(3',4'-l), 

(3', 4'-(d-l)), (3', 5'), (4'-l, 4'-d), . . . , (4'-(d-l), 4'-d), (4'-d, o), 

(5,o)}. 

Since 04 has delay d, we simply follow Lemma [T] and add d 
neurons in II corresponding to 04. Let t be the time when 0-3 
and 03/ spike. The time that o spikes the second time (since 
the spike from 0-5 makes o D spike the first time, followed by 
the delayed spike from 0-4) i.e. the halting time, is t + d + 1 
and the environment receives two spikes in total. This time 
coincides with the halting time of IT, which also sends two 
spikes to its environment. and R 2 are both satisfied for 
this case. 



In the case where both 04 and 05 have delays di and ds 
respectively, then following Lemma [1] II has di + ds addi- 
tional neurons. For both systems, halting time is t+d m ax + 1 
where d max = max(d4,ds), and two spikes are sent to the 
environment, satisfying R' x and R^. □ 

We can now have the following theorem. 

Theorem 1. Given an SNP system II with delays con- 
taining one or more of the following routings: sequential, 
iteration, join, split, there exists an SNP system II that sim- 
ulates n. 

PROOF. Proof follows from Lemma [T] [2] El Q] □ 

Notice that Observation [T] and [5] hold for all four routing 
constructs. 

4. FINAL REMARKS 

We have presented an alternative construction of a II that 
simulates a given II, improving the previous work so that 
we use only one initial spike for both systems. The halting 
time of both systems also exactly coincide with one another. 
The trade off is that there is an "explosion" of neurons in II 
for every delay in LI i.e. we add d; neurons in II for every Oi 
that has a delay. 

For our further work, we will consider nondeterministic SNP 
systems i.e. neurons having more than one applicable rule, 
since we only consider deterministic systems in this work. 
Minimization of the number of neurons of a LI simulating a II 
is also desirable, including providing bounds to the number 
of neurons, spikes, and types of spiking rules. We will also 
use the matrix representation of SNP systems without delays 
from [15] and then use massively parallel processors such as 
graphics processing units to create simulations of computa- 
tions as was done in [2]. Lastly, certain results and appli- 
cations of SNP systems that use delays can be converted to 
SNP systems without delays e.g. generating automatic se- 
quences as in [4], and performing arithmetic operations as 
in [16] . among others. 
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